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We study patterns of chiral symmetry breaking at zero temperature and its sub- 
sequent restoration at nonzero temperature within the SU{3)r x SU{3)e linear sigma 
model. Gap equations for the masses of the scalar and pseudoscalar mesons and the 
non-strange and strange quark condensates are systematically derived in the Hartree 
approximation via the Cornwall- Jackiw-Tomboulis formalism. In the chiral limit, 

■ the chiral symmetry restoring transition is found to be first order, as predicted by 
' universality arguments. Taking the experimental values for the meson masses, how- 
, ever, the transition is crossover. The absence of the U{1)a anomaly is found to drive 
' this transition closer to being first order. At large temperatures, the mixing angles 

between octet and singlet states approach ideal fiavor mixing. 

CJ ; I. INTRODUCTION 

Ch . . . 

Chiral symmetry is broken in the vacuum of quantum chromodynamics (QCD). At temperatures 

of about 150 MeV, lattice QCD calculations indicate that chiral symmetry is restored |^]. The order 
, of the phase transition depends on the mass of the non-strange up and down quarks, m„ ~ to^, and 

■ the mass of the strange quark, [^. In nature, to„ ~ rud '-^ 10 MeV and rUs ~ 100 MeV ||]. At 
temperatures on the order of 150 MeV, heavier quark flavors do not play an essential role. 

For Nf massless quark flavors, the QCD Lagrangian has a SU {Nf)r x SU {Nf)i x U{1)a symmetry. 
In the vacuum, a non- vanishing expectation value of the quark condensate, {q£ qr) ^ 0, spontaneously 
breaks this symmetry to the diagonal SU{Nf)v group of vector transformations, V = r + £. For 
Nf = 3, the effective, lovif-energy degrees of freedom of QCD are the scalar and pseudoscalar mesons. 
Since mesons are qq states, they fall in singlet and octet representations of SU(3)v- 

The SU(Nf)r x SU{Nf)e x U(1)a symmetry of the QCD Lagrangian is also explicitly broken by 
nonzero quark masses. For M < Nf degenerate quark flavors, a SU{M)v symmetry is preserved. 
If M < Nf, the mass eigenstates are mixtures of singlet and octet states. For instance, in the 
pseudoscalar meson sector, this mixing occurs between the rj and the rj' meson, with the 77 meson 
being mostly octet and the rj' meson being mostly singlet, with a mixing angle of about -10 to -20 
degrees [||. 

As shown by 't Hooft instantons also break the U(1)a symmetry explicitly to Z{Nf)A @- For 
the low-energy dynamics of QCD, however, this discrete symmetry is irrelevant. 
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Pisarski and Wilczek |g] discussed the order of the chiral transition using renormahzation group 
arguments in the framework of the hnear sigma modeL This model is the effective theory for the 
low-energy degrees of freedom of QCD and incorporates the global SU{Nf)r x SU{Nf)i x U{1)a 
symmetry, but not the local SU{3)c color symmetry. They found that, for Nf = 2 flavors of 
massless quarks, the transition can be of second order, if the U{1)a symmetry is explicitly broken 
by instantons. It is driven first order by fluctuations, if the U{1)a symmetry is restored at Tc. For 
Nf = 3 massless flavors, the transition is always first order. In this case, the term which breaks the 
U{1)a symmetry explicitly is a cubic invariant, and consequently drives the transition first order. 
In the absence of explicit U{1)a symmetry breaking, the transition is fluctuation-induced of first 
order. 

In nature, the chiral symmetry of QCD is explicitly broken by nonzero quark masses. In this 
case, one has to resort to numerical calculations to determine the order of the chiral transition. At 
present, however, lattice QCD data have not unambiguously settled this issue. For physical values 
of the quark masses, calculations with staggered fcrmions favor a smooth crossover transition, 
while calculations with Wilson fermions predict the transition to be first order. 

As an alternative to lattice QCD calculations, one can also use the linear sigma model to make 
predictions on the order of the phase transition in QCD. Furthermore, various symmetry-breaking 
scenarios can be more easily investigated than on the lattice. Studying the linear sigma model at 
nonzero temperature, however, requires many-body resummation schemes, because infrared diver- 
gences cause naive perturbation theory to break down Q . 

For Nf = 2, effective chiral models for QCD have been studied extensively, because in this case 
SU{2)r X SU{2)i is isomorphic to 0(4), and 0{N) models Q in general are particularly amenable 
to many-body approximations at nonzero temperature. For an incomplete list of references, see 
[0, where some of us (J.T.L. and D.H.R.) studied the 0{N) model in the Hartree and large-iV 
approximations. 

For Nf — 3, many-body resummation schemes become considerably more involved due to the 
larger number of degrees of freedom. The SU(S)r x SU{3)i linear sigma model was previously 
studied at nonzero temperature in fl^ . The genuine problem of these approaches is that they employ 
methods related to the standard loop expansion to compute the effective potential. At nonzero 
temperature, the loop expansion is known to fail in the case of spontaneously broken symmetry, 
since it generates imaginary masses for the particles. The physical reason for this failure is that 
contributions from thermal excitations to the masses are neglected. This can be amended by self- 
consistent resummation schemes like the mean-field or the Hartree approximation. In [|l3|, one of 
us (J.S.-B.) studied the linear sigma model in the mean- field approximation. 

In this work, we study the SU{3)r x SU{3)i linear sigma model in the Hartree approximation. 
We derive this approximation systematically from the Cornwall-Jackiw-Tomboulis (CJT) formalism 
[0. Thus, we extend previous work for Nf = 2 [|o| to Nf = 3. We study possible patterns of 
symmetry breaking in the vacuum and its subsequent restoration at nonzero temperature. We 
focus on both the cases where chiral symmetry is and where it is not explicitly broken by nonzero 
quark masses. Lattice QCD data indicate that the U{1)a anomaly becomes small near the chiral 
transition [T^ . Therefore, we also study the influence of the U{1)a anomaly on symmetry restoration 
at nonzero temperature. 

This paper is organized as follows. In Sec. ||, we introduce the SU{3)r x SU{3)e linear sigma model. 
In Sec. |l^ the possible patterns of symmetry breaking in the vacuum are discussed. Section |^ is 
devoted to the vacuum properties of the model at tree level. In Sec. 0, we derive the gap equations for 
the condensates and the masses in Hartree approximation via the CJT effective potential. Numerical 



results are presented in Sec. VI. We conclude this work in Sec. [VII| with a summary of our results 



We use the imaginary-time formalism to compute quantities at nonzero temperature. Our notation 

is 

jj{k)^T J ^f{2mnT,k) , J f{x)^ dr J d'^ f{r,^) . (1) 

We use units /i = c = fcs = 1. The metric tensor is g^" = diag(+, — , — , — ). Throughout this work, 
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all latin subscripts are adjoint U{3) indices, a = 0, ... ,8, and a summation over repeated indices is 
understood. 



II. THE LINEAR SIGMA MODEL FOR THREE FLAVORS 

The Lagrangian of the SU{3)r x SU{3)i linear sigma model is given by 

£($) = Tr - $t$) _ [Tr ($1'$)]^ - AsTr 

+ c[Det($)4-Det($t-)] +Tr[iJ($ + $t)] . (2) 

$ is a complex 3x3 matrix parametrizing the scalar and pseudoscalar meson nonets, 

^^Ta(t>a^Ta{<Ja+ina) , (3a) 

where Ca are the scalar fields and tTq are the pseudoscalar fields. The 3x3 matrix H breaks the 
symmetry explicitly and is chosen as 

H = Taha , (3b) 

where ha are nine external fields. Ta — Aa/2 are the generators of U{3), where Aa are the Gell-Mann 
matrices with Aq — ^^1- The Ta are normalized such that Tr(TaTf,) — Sabf^ and obey the C/(3) 
algebra with 

[Ta,Tb] ^ i fabcTc , (4a) 

{Ta.Tb} ^ dabcTc , (4b) 

where fabc and dabc for a, 6,c — 1, ... ,8 are the standard antisymmetric and symmetric structure 
constants of SU{3) and 

fabO = , dabO = Sab ■ (4c) 

In Eq. (||), is the tree- level mass of the fields in the absence of symmetry breaking, Ai and 
A2 are the two possible quartic coupling constants, and c is the cubic coupling constant. In four 
dimensions, the cubic and the two quartic terms are the only relevant SU{3)r x SU{3)i invariant 
operators. 

The terms in the first line of Eq. (|^) are actually invariant under the larger group of U{3)r x U{3)f, 
symmetry transformations, 

4> — yUr<PUl , Ur,i = exp {i uj^j T") . (5) 

Introducing = (a;°±w")/2, the right- and left-handed symmetry transformations can be 

alternatively written as vector, V = r + £, and axial vector, A — r — £, transformations. It is 
then obvious that $ is a singlet under U{l)v transformations exp{iujyT^). This U{l)v is the U{1) 
of baryon number conservation and thus always respected. The terms in the first line of Eq. (^ are 
therefore invariant under SU{3)r x SU{3)e x U{1)a = SU{3)v x U{3)a- 

The determinant terms correspond to the U{1)a anomaly in the QCD vacuum. As shown by 
't Hooft they arise from instantons. These terms are invariant under SU(3)r x SU{3)e ^ 
SU{3)v X SU{3)a transformations, but break the U{1)a symmetry of the Lagrangrian explicitly. 
The last term in Eq. (0) breaks the axial and possibly the SU{3)v vector symmetries explicitly. The 



patterns of explicit symmetry breaking will be discussed in detail in Sec. III. 

The (Ta fields are members of the scalar {J^ = 0^) nonet and the iTa fields are members of the 
pseudoscalar (J^ = 0~) nonet, 
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T,aa = ^| _-i=a>; + -i=a8 + ^ao ], (6a) 
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Here, tt* = (tti ± z7r2)/\/2 and tt" = tts are the charged and neutral pions, respectively. 
(7r4 ± j7r5)/\/2, = (TTfi + i 7r7)/\/2, and iiT" = (ttb - i 7r7)/\/2 are the kaons. In general, because 
the strange quark is much heavier than the up or down quarks, the ttq and the TTg are admixtures of 
the r] and the 77' meson. 

The situation with the scalar nonet is not as clear. The parity partner of the pion is the ao(980) 
meson, i.e., = [ax ± i u-i)! \fi and aj] = u^. We identify the parity partner of the kaon with the k 
meson [now referred to as Xo(1430) in [Q]. Finally, in general the ctq f^nd the a% are admixtures of 
the a [now also referred to as /o(400— 1200)] and /o(1370) mesons. [Instead of the /o(1370) meson, 
one could have chosen the /o(980) meson. However, as we shall see in Sec. IV, in the linear sigma 
model the mass of this state is closer to the /o(1370).] For details concerning the phenomenological 
status of the scalar nonet, see, for instance [l^ . 

In principle, there is the possibility that the scalar particles are not diquark states, but formed 
from two quarks and two antiquarks [ p^ . Then we would associate the k, with the k(900) discovered 
in 'kK scattering [ p8| . Determining the vacuum properties of the linear sigma model, cf. Sec. 
the mass of the k turns out to be between that of the k(900) and the iiro(1430), while its width |^] 
is closer to the observed width of the k(900). 

Symmetry breaking gives the $ field a vacuum expectation value, 

($) =Taaa. (7) 
Shifting the $ field by this vacuum expectation value, the Lagrangian can be rewritten as pO| ], 

4 \ ( ^ \ 

Gabc ~ -^J^abcdl^d j (TaC^bf^c—'i [Q abc + T^'Habcd (^d \ T^aT^bf^c 
— 2Habcd<7aO'bTrcT^d — -^^abcdiCaCb^cC^d + T^aT^bT^cT^d) - U{a) , (8) 

where 

m,'^ 1 

U{a) — — ~ Gabc^a^b^c + -^^abcd^a^b^c&d — ha^a (9) 

is the tree-level potential and da is determined on the tree level by 

dU{a) 2 - 4^ , „ 

^ _ = TO CTa - 3 yabcCTbCTc + T:-^ abcdCrbUcCrd - ha ^ . (10) 
0(7a 3 

The coefficients Gabc, ^abcd, and Habcd are given by 



Gabc — ^ 
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dabc — 2 i^aodobc + SbodaOc + ^codabo) + 2 '^OOO'^aO'^bO^cO 



(11a) 



^abcd = {Sabred + ^ad^bc + ^ac^bd) + [dabndncd + dadndnbc + dacnd-nbd) , (Hb) 

Habcd — —r 6ab5cd H~ Q {dahndjicd H~ facnfnbd fbcnfnad^ ■ (-^^^) 
4 o 
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The tree-level masses, {m'g)ab and (mp)afc are given by 

{ml)ab = Sab - 6 Gabc ^0 + -^ J^abcd ^c^d , (12a) 

{my)ab Sab + 6 Qabc + 4 Habcd ^c^d ■ (12b) 

In general, these mass matrices are not diagonal. Consequently, the fields (ctq, tTq) in the standard 
basis of C/(3) generators are not mass eigenstates. Since the mass matrices are symmetric and real, 
diagonalization is achieved by an orthogonal transformation, 

= ^ , (13a) 

7^, = O^^^ TTa , (13b) 



III. PATTERNS OF SYMMETRY BREAKING 



In this section, we discuss possible patterns of symmetry breaking in the vacuum. We begin 
with the most symmetric case, i.e., with the minimum number of nontrivial couplings, and then 
successively reduce the symmetry. 

1. if = 0, c = 0, A2 = 0: For > 0, the symmetry group is 0(18), on account of 

Tr($t$) = i(a^-f^^) . (14) 

The physics of the 0{N) model has been studied extensively in the past and so we shall 

restrict ourselves in the following to A2 0. Here, we only mention that for < 0, the 0{18) 
symmetry is spontaneously broken to 0{17) and there are 17 Goldstone bosons. 

2. H — 0, c — 0, X2 0: For rri^ > 0, the Lagrangian has a global SU{3)v x U{3)a symmetry. 
For < 0, Eq. (^ shows that <I> develops a non- vanishing vacuum expectation value, and the 
U{3)a symmetry is spontaneously broken. By the Vafa-Witten theorem |2l| , only the axial 
symmetries can be spontaneously broken, while the vector symmetries remain intact. One can 
distinguish two cases |Q: 

(a) A2 > 0. SU (3)y X U{3)a is broken to SU{3)v, with ($) ^ diag(l, 1,1) and the appearance 
of 9 Goldstone bosons which comprise the entire pseudoscalar nonet, i.e., the pions, the 
kaons, the rj, and the 77' meson. The nine massive scalar particles fall into irreducible 
representations of SU{3)v- Since the mesons consist of a quark (a [3] of SU{3)v) and 
an antiquark (a [3]), these representations are a singlet and an octet, because [3] x [3] = 
[1] -t- [8]. The mass of the singlet, the a meson, is in general different from the (degenerate) 
masses of the octet particles. 

(b) A2 < 0. SU{3)v X U{3)a is broken to SU{2)v x U{2)a, with ($) - diag(0,0, 1) and 10 
Goldstone bosons. 

3. H — 0, c 0, X2 0: The symmetry is SU{3)v x SU{3)a- A non- vanishing ($) spontaneously 
breaks this symmetry to SU{3)v, with the appearance of 8 Goldstone bosons, which is the 
complete pseudoscalar octet. The ninth Goldstone boson of case 2(a), the rj' meson, becomes 
massive and thus is no longer a Goldstone boson, because the U{1)a symmetry is already 
explicitly broken. The masses of the scalar particles behave as in case 2(a). Note that from 
Eq. (0), < is no longer required for spontaneous symmetry breaking when c 7^ 0. 
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4. iJ ^ 0, c = 0, A2 ^ 0: In QCD, this corresponds to non-vanishing quark masses, but vanishing 
U{1)a anomaly. Since ($) must carry the quantum numbers of the vacuum, only hg, /13, and 
hs can be nonzero. One can distinguish three cases: 

(a) ho =/= 0, /i3 = hg = 0. All quark masses are equal, i.e., m„ = = m^. In this case, 
the U{3)a axial symmetry is explicitly broken, i.e., the 9 Goldstone bosons of case 2(a) 
become (mass degenerate) pseudo- Goldstone bosons. The SU(3)v symmetry remains 
intact and the scalars follow the classification as discussed in case 2(a). 

(b) ho 7^ 0, /13 = 0, hs 7^ 0. Only the non-strange flavors are degenerate in mass, i.e., 
ruu = rrid ^ rUs- Besides the explicitly broken U{i)A symmetry, SU{'i)v is explicitly 
broken to SU{2)v- For the scalar particles, the following applies. If there was ideal 
flavor mixing, i.e., the physical particles are also eigenstates of flavor, one particle is 
an ss state (the /o meson), while all others contain at least one non-strange quark or 
antiquark. The latter then fall into irreducible representations of SU{2)v- For the scalar 
particles containing no strange quark or antiquark, a quark (a [2] of SU{2)v) couples 
with an antiquark (a [2]) to form the tr meson singlet and the ao meson triplet, since 
[2] X [2] = [1] -I- [3]. The strange scalar particles have only one quark in a [2] or a [2] 
representation, and therefore fall into doublets of SU{2)v- The k"*" and k° mesons form a 
[2], while the k~ and mesons form a [2]. Because the masses of quarks and antiquarks 
are identical, these two doublets are mass degenerate. In nature, however, flavor mixing is 
not ideal and the /o meson has a qq admixture, just as the a meson has an ss admixture. 
For the pseudoscalars, only the four non-strange pseudo-Goldstone bosons, the pions and 
the r]' meson, are degenerate in mass. The kaons and the r] meson both have different 
masses. Since the pions are non-strange qq states, and the 77' meson is degenerate in mass 
with the pions, it follows that it is also a non-strange qq state. Then, the 77 meson is a 
pure ss state, i.e., flavor mixing is ideal in the pseudoscalar sector. 

(c) ho 7^ 0, /13 7^ 0, /i8 ^ 0. Here, SU{2))v is completely broken. Even the non-strange 
pseudo-Goldstone bosons are no longer completely degenerate in mass. 

5. g 7^ 0, c 7^ 0, A2 7^ 0: Now, from the f/(3)A = SU{i)A x U{\)a symmetry, the U{1)a is ex- 
plicitly broken by instantons. Again, there are three cases: 

(a) ho ^ 0, h^ — hg — 0. In this case, the remaining SU{i)A axial symmetry is explicitly 
broken, i.e., the 8 Goldstone bosons of case 3 become (mass degenerate) pseudo-Goldstone 
bosons. As above, the SU{2>)v symmetry remains intact. The scalar particles behave as 
in case 4(a). 

(b) ho 7^ 0, /13 = 0, /i8 7^ 0. Besides the explicitly broken SU{2>)a symmetry, SU{2>)v is 
explicitly broken to SU{2)v- The scalar and pseudoscalar particles behave as in case 
4(b), except that the ry' meson mass is different from the pion mass because of the U{1)a 
anomaly. Then, flavor mixing is no longer ideal in the pseudoscalar sector. 

(c) ho ^ Q, h^i, ^ 0, hfi 7^ 0. This is the case realized in nature, although violation of isospin 
SU{2)v is small (the charged and the neutral pions are almost degenerate in mass). 

In the following, we shall restrict ourselves to studying cases 2(a), 3, 4(b), and 5(b). The first 
two cases are interesting because they represent the idealized scenario where the quark masses are 
zero, i.e., the chiral limit. The last two cases are close to the situation in nature where quark masses 
break the chiral symmetry explicitly. Lattice QGD data indicate that the U{1)a anomaly becomes 
small for large temperatures |15). This motivates our interest in cases 2(a) and 4(b), see also [ p^ . 
Since isospin SU{2)v violation is rather small in nature, it is sufficient and easier to study case 5(b) 
instead of 5(c). 
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IV. CONDENSATES AND MASSES IN THE VACUUM 



In this section, we establish the vacuum properties of the SU{3)r x SU{3)e model in the various 
cases selected above. In contrast to the previous section, it is now easier to begin with the most 
asymmetric case 5(b). This is the case closest to nature and therefore it is natural to use the 
experimental values for the meson masses and decay constants as input to determine the coupling 
constants. The linear sigma model has six parameters, m^, Ai, A2, c, ho, and hg. It therefore 
requires six experimentally known quantities as input. We choose to^, ttik, J-k, Ik, the average 
squared mass of the ?/ and ry' mesons, + m^;, and m„. The other masses, mag, m^, m/^, the 
difference — m^/, and the mixing angles are then predicted. 

In the other cases, nature does not provide us with information about the masses and the decay 
constants. In case 4(b), c = 0, and we need to specify only five input parameters. It turns out that 
we do not need to fix + m^, , because m^/ = TOtt , and the mass of the 77 meson is predicted as 
well. 

The cases 2(a) and 3 correspond to the chiral limit which is not realized in nature. Therefore, 
in principle one should not use experimental values for the masses and decay constants to fix the 
parameters m^, Ai, A2, and, in case 3, c. For the sake of dcfiniteness, however, we use the following 
tentative generalizations of experimental data. We use the pion decay constant, /tt, extrapolated to 
the chiral limit, the a meson mass, and, in case 3, the t]' meson mass. [In case 2(a), m^/ — = 0.] 
For the final input parameter, since the scalar octet is degenerate in mass when the SU (3)^ x SU (3)^ 
symmetry is not explicitly broken, we use an average of the experimental values for the masses of 
the scalar octet. 



A. Explicit chiral symmetry breaking with U{1)a anomaly 

The vacuum expectation value is ($) = Tq o'o + Js The equations for the condensates (Tq and 
as read 



ho = 
hs = 



- CTo + ( Ai + ^ ) CT^ 



Co 



2%/6 
A2 



+ (Ai + A2) (To 



™ 71 271 ^ + A2) CTq - O-O CTs +( Ai + ^ ) cr. 



A2 _ 




^3^2"*^ 













(15a) 
(15b) 



The PCAC relations (see Appendix) determine the values of the condensates from the pion and kaon 
decay constants, /^r, /k, 



Co 



V6 



o's = (/tt - Ik) ■ 

We use the experimental values f-^ = 92.4 MeV, Jk — 113 MeV 
The nonzero elements of the scalar mass matrix are 



(16a) 
(16b) 



(™l)oo = 



I cCTo + (3Ai + A2) ffg + (Ai + A2) al 



(™s)ii = ("^5)22 = ('"s) 



33 



(To ^ 0-8 + (Ai + A2)o-o + \/2 AsCToCTs + (^1 + y) '^s > 



V6 

{^Dm = ("^5)55 = ("^s)66 = 



77 



(17a) 
(17b) 
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A 



0-8 + (Ai + A2) CTq —ao(Ts+[Xi + —]a\ 



V2 



V3 

-^ + 2(Ai + A,)ao-^a8 



(mDgg = + (To + -7^ cTs + (Ai + \2)<7l - V2X20-oas + 3 ( Ai 
(™l)o8 = i^^s) so = 



A, 



-2 

^8 ' 



0-8 



(17c) 
(17d) 
(17e) 



While the masses of the qq and the k mesons are given by the (11) and (44) elements of the mass 
matrix, = (?ti|)ii, to^ = (7715)44, the a and /o meson masses are obtained by diagonalizing the 
(08) sector of the mass matrix. According to Eq. (p^, 

ml EE (m|)o = (m|)oo cos^ 6*5 + (m|)88 sin^ 6*5 + 2 (to|)o8 cos 6*5 sin 6*5 , 
m? = (m|)8 = (m|)oo sin^ 0s + {ml)s8 cos^ - 2 (to|)o8 cos 6*5 sin 6*5 , 



where the scalar mixing angle 6s is given by 

tan 26*5 



(18a) 
(18b) 



2(m|) 



5^08 



("i^)oo - (mi;) 



88 



The pseudoscalar mass matrix is given by 



'PJii 



'■P)22 



'33 



V2 



do + 0-8 +( Ai + — ) CTq + — ^2 CTo 0-8 +( + — ) CT, 



'PJ66 

c 



'77 



V6 



2V3 

c 



*8 + Ai + 



CTo + Ai + - A2 0-, 



CT8 + Ai + 



'80 



c 2 ^ 



3%/2 

— A2 (To 0-8 + Al + 



6 
A2 



3V2 



0-8 



(18c) 

(19a) 
(19b) 

(19c) 
(19d) 
(19e) 



While the mass of the pion and the kaon are given by the (11) and (44) elements of the mass matrix, 
toJ = (mp)ii, = (mp)44, the rj' and r/ meson masses are obtained by diagonalizing the (08) 



sector of the mass matrix. According to Eq. (y_3|), 

mfy = (mp)o = {iTip)oo cos^ 9p + {mp)s8 sin^ 9p + 2 (mp)o8 cos9p sin^p 
mf^ = (mp)8 = (TOp)oo sin^ 9p + {m%)ss cos^ 9p — 2 (mp)o8 cosf?p sinflp 

where the pseudoscalar mixing angle 9p is given by 



tan2^?p — 



2{mj,) 



08 



(20a) 
(20b) 

(20c) 



(™p)00 - (™p)88 

Th e ex plicit symmetry breaking terms, ho and hs, are determined from Eqs. ([Tsl), (p^), ( |l9b| ), 
and (19c) as 



^0 = {ml U + 2ml^ fx) , 
2 

^8 = {ml U - m\ fx) ■ 



(21a) 
(21b) 
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Using the experimental pion and kaon masses, one obtains Hq = (286.094 Me V)'^ and hg = 
-(310.960 Me V)3. 

Comparing Eqs. (15b) and (17c) reveals 

hs = ml<Js , (22) 

i.e., the mass of the k meson is predicted to be m„ = 1124.315 MeV, which is about 21% smaller 
than the experimental value. The average ry and rj' meson mass squared determines A2 

_ 3(2/k-A) ml-{2fK + U) m^-2 (mg,+m^) (/k - ^) 

[3f^ + 8fK{fK-.U)]{fK~f.) ■ ^ ' 

For m^, + TO^ = (1103.625 MeV)2, one obtains A2 = 46.484. 

The difference of the pion and kaon masses squared and A2 determine c. 



Ik - /tt 

Now, also the mass of the ao meson is fixed. 



A2 (2 fx - U) = 4807.835 MeV . (24) 



m 



^rnl + [Jk - U) [c - A2 (2 + /„)] = (1028.707 MeV)^ , (25) 



which is about 4% larger than the experimental value. The pseudoscalar mixing angle is 9p — —5°. 
This angle determines the individual rj and 77' meson masses to be m,, — 539.008 MeV and to,,/ = 
963.046 MeV. These values are surprisingly close to the experimental ones; the rj meson is about 2% 
lighter and the rj' meson is about 0.6% heavier than in nature. 

Finally, Ai is determined by fixing either the mass of the a or the /o meson; the other mass, the 
scalar mixing angle, 63, and are then given by solving a nonlinear equation for the fixed mass. 
Here, we choose rria = 600 MeV to yield Ai = 1.400 and to/q = 1221.113 MeV, about 11% smaller 
than the experimental value. The scalar mixing angle is 9s — 19.859°, and ~ (342.523 Me V)^. 



B. Explicit chiral symmetry breaking without U{1)a anomaly 



In the absence of the U{1)a anomaly, c = 0, the condensate equations (15a) and (15b) and the 
equations for the scalar and pseudoscalar mass matrices (|l^) and (jl^) simplify. Equations (|l6| ) and 
(pT|), however, remain the same and thus yield the same values for ctq, cts, ho, and hg as above. The 
mass of the k meson is still given by Eq. ( p^ ) and is the same as above. 

The differences begin with the mass of the rj' meson, which is now identical to the pion mass, 
since the U(1)a anomaly is absent. Furthermore, A2 is given by the kaon and pion masses and decay 
constants, 

This then determines the mass of the rj meson, 

ml = ml + 2X2 !k {Ik -.!.) = (634.818 MeV)^ , (27) 

and the ao meson, 

ml, = ml - A2 (2 fl -UfK- fl) = (850.387 MeV)^ . (28) 

The pseudoscalar mixing angle is given by tan26'p = 2\/2, or Op = 35.264°. This corresponds to 
ideal flavor mixing. 

As before, Ai is given by solving the equation for the a meson mass, which in turn yields the 
/o meson mass, the scalar mixing angle, and m^. For m„ = 600 MeV we find Ai — —4.550, 
m/o = 1341.367 MeV, 63 = 31.326°, and = -(503.551 Me V^). 
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C. Chiral limit with U{1)a anomaly 



In this case, ($) = To (Tq. Using fe g — and the exphcit form for the Qabc and Tabcd, the equation 
for the condensate at the tree level (15a) simplifies to 



= 



m - -j=ao + y\i + — ] (Jq 



(29) 



From the PC AC relations (see Appendix), we now obtain (Tq = ^/3/2 Lattice QCD data indicate 
a linear behavior of with the quark mass, /^r — arrig + h [E3|. Extrapolating these data to the 
chiral limit, mq (as appropriate for H — 0), one obtains /^r = 90 MeV. 

The scalar and the pseudoscalar mass matrices are diagonal and have the particularly simple form 



TO^ + ^ CTq + (Al + A2) O-Q 

v6 



5ab 



V6 



^0 + Al + 



A, 



2 \ ^2 




- c — 2Aio-o 1 ctq (5ao5fco 



CfToSaoSbQ . 



(30a) 
(30b) 



On account of Eq. (|29|), the pseudoscalar octet, a, 6 = 1, . . . , 8, is massless, as expected from group 
theory, see case 3 in Sec. III. The singlet, the 77' meson, is massive, because the U{1)a symmetry is 
explicitly broken by the U{1)a anomaly. 



,2 _ 



If we use the experimental value of the rj' meson mass, c — 6791.157 MeV. 
The scalar particles fall into a singlet, the a meson, with mass 



(31) 



4)c 



= m2-cA + -(3Ai + A2)/^ , 



(32) 



and an octet comprising the ao, k, and fa mesons. 



,2 _ 



(33) 



where the last identity follows from Eq. ([291). This equation determines Ai and A2 for given rricr, 
rUaQ, and m^'. As the octet is mass degenerate, we average the experimental values for "i^^, m^, and 
m^^, weighted by the respective isospin degeneracy, to obtain an average mass of ma„ = 1225.795 
MeV. For the a meson we again take = 600 MeV. This results in Ai — —9.291 and A2 = 110.046. 
Then, follows from any of the mass equations as = —(164.921 MeV)^. 

Note that the scalar singlet is lighter than the scalar octet, but the pseudoscalar singlet is heavier 
than the pseudoscalar octet. This "inverted mass spectrum" for the scalar mesons relative to the 
pseudoscalar mesons is a general feature in the presence of the U{1)a anomaly. It arises from 
the relative difference in sign of the terms ~ Qabc in Eqs. (f2|). 



D. Chiral limit without U{1)a anomaly 



In this the U{1)a symmetry is not explicitly broken, all nine pseudoscalar particles are 

massless. This is readily derived from Eqs. (29) and ( |301:| ) when c = 0. 
The scalar particles again fall into a singlet, the a meson, with mass 



K) 



00 



m2 + -(3Ai+A2)/2 



(34) 
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and an octet comprising the ag, k, and /o mesons, 

^lo ^ i^Dii - m2 + ^ (Ai + A2) fl ^ ml - iX^fl ^X2fl. (35) 

where the last identity follows from Eq. ( p9| ) with c = 0. This equation determines Ai and A2 for given 
rria and TOa,,. We again use = 1225.795 MeV and m<^ = 600 MeV. This results in Ai = -47.019 
and A2 = 185.503. Then, follows from any of the mass equations, m? = — (424. 264 Me V)^, which 
is negative, as required in this case for spontaneous symmetry breaking. 



V. THE EFFECTIVE POTENTIAL IN THE CORNWALL-JACKIW-TOMBOULIS 

FORMALISM 

The effective potential of the Cornwall-Jackiw-Tomboulis formalism |l^ is 



2 Jk 

Here, U{a) is the tree-level potential of Eq. (||), and 

Sabik;ij)^~k'Sab+{ml)^^ , (37a) 
Pab'ik; a) = -e Sab + {ml)^^ , (37b) 

are the tree-level propagators for scalar and pseudoscalar particles, with the respective mass matrices 
®- 

The expectation values for the scalar fields, Ua, and the full propagators for scalar, S{k), and 
pseudoscalar, V{k), particles are determined from the stationarity conditions 

^0, (38a) 



SSab{k) 

SV[<j,S,V] 



5Vab{k) 

With Eq. (p6|), the latter two can be written in the form 



, (38b) 
. (38c) 



where 



S:b\k) = S~,\k-a) + Y.,b{k) , (39a) 
Kb ik) = Pab' (fc; ^) + ^abik) , (39b) 



(40a) 



SSbaik) 
SVbaik) 



n..(fc)- 2^^114^, (40b) 
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are the self-energies for the scalar and pseudoscalar particles. 




FIG. 1. The double-bubble diagrams. Full lines are scalar particles, dashed lines are pseudoscalar particles. 

In general, V2[d',S,'P] is the sum of all two-particle irreducible (2PI) diagrams, with all lines 
representing full propagators. Here, we restrict ourselves to the most simple class of 2PI diagrams, 
the double-bubble diagrams of Fig. |l], which is equivalent to the Hartree approximation. Explicitly, 



V2[a,S,'P] = Tabcd 



Sabik) / Scdip) + / Vabik) / Vcdip) 
k <J p 'J k <J p 



+ 2nabcd j^Sab{k) J Vadip) • (41) 



Note that, in the Hartree approximation, V2 is actually independent of aa- Therefore, the stationarity 
conditions for the condensates are 



ha =771^ Ua- iQabc jcTbCTc + j [Scb{k) - Vcbik)] 



+ '^^abcd 



■ CTb 



Scb{k) (Jd + 'iHbcadCTd jPcbik) 
Jk J Jk 



(42a) 



Since the self-energies ( ^0| ) are independent of momentum in the Hartree approximation, the 
Schwinger-Dyson equations (p9|) for the full propagators assume the simple form 



Sabik) = -kHab+{Ml)^^ , (42b) 
Kbik) = -kHab+{Ml)^^ , (42c) 

where the scalar and pseudoscalar mass matrices in the Hartree approximation are given by 



mab 



Sab 


— 6 Qabc Cc - 


1- 4 !Fabcd 


fc <Jd + 


/ Scdik) 

Jk 


+ 4 Habcd / 
Jk 


rr? Sab 


+ 6 Qabc CTc - 


\- 4 Habcd 


CTc ^d 4 


- j^S,d{k) 


^^Tabcd f 
Jk 



Vcd{k) 



(43a) 
(43b) 



In general, the mass matrices are not diagonal in the standard basis of U{i) generators, see Sec. IV. 
Consequently, the propagators are also not diagonal in this basis. Physically, however, only mass 
eigenstates can propagate. Th erefore, we have to diagonalize the propagators before we compute 
the loop integrals in Eqs. (42a) and (^). 

In the Hartree approximation, all particles are stable qu asip article s, i.e ., the imaginary parts of the 
self-energies vanish. Therefore, the inverse propagators (42b) and (42c) are real- valued. They are 
also symmetric in the standard basis of f7(3) generat ors a nd thus diagonalizable via an orthogonal 
transformation. This transformation is given by Eq. dT3q), with the obvious replacements 



(44) 
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In cases 2(a) and 3 in Sec. 11], the mass matrices are diagonal at zero temperature, cf. Sec. IV, 
and are taken to be diagonal at nonzero temperature T as well. In cases 4(b) and 5(b), they have 
ofF-diagonal elements in the (08) sector at zero temperature, and consequently also have ofF-diagonal 
contributions in this sector at nonzero temperature. In the latter case, diagonalization proceeds as 



in Eq. (13c) with the replacements ([44|) 



The propagator matrices are diagonalized by the same orthogonal transformation as their inverse. 



The loop integrals in Eqs. ( [42a| ) and (43) are therefore computed, for example, as 

Sam = oif / sm oif^ , (45) 



fc 

where Si{k) is the scalar propagator in the mass eigenbasis. 

The loop integral in Eq. ( |45|) requires renormalization. Renormalization of many-body approxi- 
mation schemes is nontrivial Tin] , but does not change the results qualitatively. We therefore simply 
omit the vacuum contributions to the loop integrals and set 

Sm-I^, ^- fexp|M^|-lV' , (46) 



and similarly for the pseudoscalar loop integrals. Here, ek[(A^|)j] = y + {Mg)i is the relativistic 
energy of the ith scalar quasiparticle with momentum k. 

VI. RESULTS 

In this section, we discuss the numerical results at nonzero temperature for the four cases of 
interest. 



A. Chiral limit 

In Fig. ||(a), the masses are shown in the chiral limit with the U{1)a anomaly, c ^ Q. This 



corresponds to case 3 of Sec. Ill, for which the zero-temperature properties were discussed in Sec. 
IV C . Accordingly, there are eight Goldstone bosons, the three pions, the four kaons and the ry 
meson, while the rj' meson has a large mass due to the U{1)a anomaly. The scalar octet, comprising 
the three oq mesons, the four k mesons, and the /q meson, is mass degenerate, while the singlet, the 
a meson, has a different mass. 

As the temperature increases, the scalar masses decrease while the pseudoscalar masses increase. 
The mass of the Goldstone bosons (the pseudoscalar octet]increases because the Hartree approxima- 
tion does not respect Goldstone's theorem at nonzero T [Q. We do not consider this to be a serious 
shortcoming: as shown in supplementing the Hartree approximation with the Random Phase 
Approximation (RPA) cures this problem. The analogous treatment within the CJT formalism will 
be deferred to a forthcoming publication |^ . 

At some critical temperature, Tc ^ 170 MeV, there is a first order phase transition between 
the low-temperature phase where chiral symmetry is broken and the high-temperature phase where 
chiral symmetry is restored and all meson masses are equal. The exact numerical value of Tc can be 
determined by computing the effective potential. 

At the point where the condensates vanish, all masses become degenerate. The reason is that at 
this point the condensate equation enforces the constraint 

= ^ |2 [5oo(fc) - Poo(fc)] - E / - ' (47) 
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which is fulfilled when all masses are equal. 

In Fig. 1(b), the non-strange and strange condensates, tpns and (ps, respectively, are shown as a 
function of temperature. In the standard basis of U{3) generators, these are defined as 




(48) 



When ($) = To (To, V^ns = ctq, and ips = cro/VS, i.e., (fs = ifns/V^, as borne out by Fig. |(b). 

In Figs. ||(c,d), the masses and condensates are shown for the case without explicit U{1)a sym- 
metry breaking, c = . Thi s corresponds to case 2(a) of Sec. |lll[ . The zero-temperature properties 
were discussed in Sec. EVD. Now, there are nine Goldstone bosons, the three pions, the four kaons. 



the 77, and the ij' meson. The behavior of the scalar octet is similar to the previous case. 

As the temperature increases, the behavior of the scalar and pseudoscalar masses is quite similar 
to the masses in Fig. ||(a), i.e., the scalar masses decrease and the pseudoscalar masses increase, 
until they become degenerate in a first order phase transition. The critical temperature for this 
transition, however, appears to be slightly lower than for c 7^ 0. A notable difference between 
Figs. |^(a) and (c) is that the masses do not become degenerate continuously as the condensates 
vanish. The reason is that, for c = 0, the above constraint equation (^^ is absent. However, these 
phenomena are irrelevant for the thermodynamic properties of the model, as they occur in a region 
where the solutions are thermodynamically unstable and the stable solution has to be found from 
Maxwell's construction for first order phase transitions. 



with anomaly 



without anomaly 




100 200 300 

T [MeV] 



T [MeV] 



FIG. 2. The meson masses and the condensates as a function of temperature for (a,b) the case with 
U{1)a anomaly (c 7^ 0) and (c,d) the case without U{1)a anomaly (c = 0), for H = 0. In (a,c), the full 
lines are the masses of the ao, k, and fo mesons, the dotted lines represent the a meson mass, the dashed 
lines are tt, K, and 77 meson masses, and the dot-dashed lines are the ri' meson mass. In (b,d), the full lines 
represent the non-strange condensate, and the dotted lines represent the strange condensate. 
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Another interesting feature is that the mass of the a meson is proportional to the condensates. 
The reason for this is that the condensate equation and the equation for the <t mass can be combined 
to give 

Ml=^-Tom^al. (49) 



B. Explicit chiral symmetry breaking with U{1)a anomaly 

In Fig. H, we show the masses for the scalars (a), the pseudoscalars (b), the condensates (c), and the 
mixing angles (d) for explicit chiral symmetry breaking, incl uding the U{1)a anomaly. The masses 



behave according to the discussion of case 5(b) in Sec. |II|. As the temperature increases, chiral 
symmetry is restored in a crossover transition and all masses become approximately degenerate. 
The temperature range of the crossover transition is ~ 220 MeV, i.e., about 50 MeV higher than 
in the case where there is no explicit chiral symmetry breaking. A notable feature is that the k 
becomes lighter than the ao meson, and the rj' meson lighter than the kaon, at about 240 MeV. 

The non-strange and the strange condensates, ip^s and ips, are shown in Fig. ||(c). Note that the 
strange condensate decreases much more slowly with temperature than the non-strange condensate. 
Figure ^(d) shows the scalar and pseudoscalar mixing angles. At large temperatures, these approach 
arcsin(l/v%) ~ 35.264°. From Eq. (|l|), 

[2 1 

o- = CTo = y - (To + cr8 , (50a) 
1 [2 

/o EE CTg = -— (To + Y 2 ^8 • (50b) 

On the other hand, 

Vns = CTo + CT8 , (51a) 

V's = - y^CTs ■ (51b) 

This shows that the a meson becomes an excitation of the non-strange condensate, i.e., a purely 
non-strange qq state. On the other hand, the /o meson is an excitation of the strange condensate 
and a pure ss state. Similarly, the rj' meson is purely non-strange and the rj meson is purely strange. 
This is what was referred to as ideal flavor mixing earlier. 

The transition occurs at temperatures which are not significantly larger than the strange quark 
mass. Therefore, the explicit SU{i)r x SU (3)^ symmetry breaking by the strange quark mass cannot 
be neglected and, at first, only the (approximate) SU{2)r x SU{2)e symmetry is restored. This means 
that the pion becomes degenerate with the a meson, and the oq becomes degenerate with the r]' 
meson. (The rj' is purely non-strange due to ideal flavor mixing.) However, due to the fact that the 
strange condensate decreases rather slowly with temperature, the explicit U{1)a symmetry breaking 
term ^ Gabc'^c in Eqs. ( p3| ) is not small. This causes the pion/tr meson mass still to be different 
from the ao/rj' meson mass. Only when both condensates approach zero, the U{1)a symmetry is 
effectively restored and the masses of all non-strange particles become degenerate. 

When the temperature becomes significantly larger than rus, the (approximate) SU{3)r x SU{3)i 
symmetry is restored. Then, all scalar octet states become degenerate, likewise all pseudoscalar 
octet states become degenerate. If this happened when the explicit U{1)a breaking term was still 
large, then the complete pseudoscalar octet would become degenerate in mass with the scalar singlet, 
and the scalar octet degenerate in mass with the pseudoscalar singlet. As it turns out, however, the 
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explicit U{1)a symmetry breaking becomes small around the same point where the (approximate) 
SU{3)r X SU{3)e symmetry is restored. 




FIG. 3. The meson masses, the condensates, and the mixing angles as a function of temperature for 
ho, hs ^ 0, — 0, and c 7^ 0. (a) For the scalar mesons, the full line is the ao, the dotted line is the k, 
the dashed line is the a, and the dot-dashed line is the /o meson mass, (b) For the pseudoscalar mesons, 
the full line is the pion, the dotted line is the kaon, the dashed line is the rj' , and the dot-dashed line is the 
ri meson mass, (c) The non-strange (full) and strange (dotted) condensates, (d) The scalar (full) and the 
pseudoscalar (dotted) mixing angles. 



C. Explicit chiral symmetry breaking without (7(1) a anomaly 

In Fig. ^, we show the masses for the scalars (a), the pseudoscalars (b), the condensates (c), and 
the mixing angles (d) for explicit chiral symmetry breaking in the absence of the U{1)a anomaly. 
The masses behave according to the discussion of case 4(b) in Sec. As the temperature increases, 
the chiral symmetry restoration crossover transition is much more rapid than in the previous case, 
and occurs at a slightly smaller temperature, ^ 200 MeV. A notable feature is the inverse mass 
ordering of the 1] meson and the kaon. At small temperatures and above the transition, the masses 
of the pion and the rj' meson are the same. In the temperature range from about 50 to 210 MeV, 
however, they differ. We perceive this to be an artefact of the violation of Goldstone's theorem in 
the Hartree approximation, cf. Fig. ^(c). 

The melting of the condensates. Fig. ^c), is similar to the previous case. Fig. ||(c). The mixing 
angles. Fig. ^(d), again approach ideal flavor mixing at large temperatures. The difference here , 
however, is that the ij and rj' mesons are also ideally flavor-mixed at zero temperature, cf. Sec. IV B| . 
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Due to the absence of the U{1)a anomaly, once the (approximate) SU{2)r x SU{2)e symmetry 
is restored, the pion, the ry', the oq, and the a mesons simuhaneously become degenerate in mass. 
(The f]' meson belongs to this class of non-strange particles due to ideal flavor mixing.) Once the 
temperature becomes large compared to the strange quark mass, the masses of the strange mesons 
converge with those of the non-strange mesons. 
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FIG. 4. As in Fig. Bl but for c = 0. 



VII. CONCLUSIONS 



In this work, we computed properties of the SU{3)r x 5*17(3)^ linear sigma model in the Hartree 
approximation at nonzero temperature. We first classified possible patterns of symmetry breaking, 
with special attention to the cases where the U{1)a anomaly is either absent or present, and the 
cases of zero, degenerate or nonzero, non-degenerate quark masses. We then determined the coupling 
constants from the vacuum values for the masses and the decay constants in the various cases of 
interest. We systematically derived the Hartree approximation within the CJT formalism. Within 
this approximation, we computed the masses of scalar and pseudoscalar particles, the non-strange 
and strange condensates, and the scalar and pseudoscalar mixing angles as a function of temperature. 
We checked that our results are consistent with the mean-field approximation employed in to 
compute these quantities. 

For the SU{Nf)r x SU{Nf)( model, in the case where the quark masses are zero, universality 
arguments predict the chiral symmetry restoring transition to be first order for Nj = 3 and Nf = 2 
in the absence of the U{1)a anomaly, and second order for Nf = 2 in the presence of the U{1)a 
anomaly fl. We find that the Hartree approximation correctly gives a first order transition in the 
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case Nf — 3. This is not necessarily an indication for the vahdity of this approximation, because 
earher work has shown that it incorrectly predicts a first order transition when Nf = 2 and the 
U{1)a anomaly is present The transition temperature is on the order of 170 MeV. 

As expected, when the U{1)a anomaly is absent, the rj' meson becomes a Goldstone boson for 
zero quark masses. A surprising result is that then the a meson mass is directly proportional to the 
condensate. 

For nonzero quark masses, ttIu = md ^ rris, we find the transition to be a crossover transition, 
but for vanishing U{1)a anomaly the crossover region is much more narrow than in the presence of 
the U{1)a anomaly. In the chirally symmetric phase, the mixing angles approach the situation of 
ideal flavor mixing, i.e., the a and r/ mesons are pure non-strange qq states, while the /o and the 77 
mesons are pure ss states. 

As an outlook, the present framework can be used as an alternative to lattice QCD studies Q 
to study the order of the chiral symmetry restoring transition as a function of the strange and 
non-strange quark masses [Q. Moreover, other meson properties such as the decay widths and the 
spectral functions can be self-consistently computed at nonzero temperature These proper- 

ties can be experimentally investigated in relativistic nuclear collisions, for instance at Brookhaven 
National Laboratory's Relativistic Heavy-Ion Collider. 
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APPENDIX: DERIVATION OF EQS. @ 

The infinitesimal form of the SU{3)r x SU{3)i x U{1)a symmetry transformation (^) is 

T" (f)" — > 0'^ - i [T", T*] (/)'' + iuj'X {T", T''} 4)'' . (Al) 
For axial- vector transformations, uiy = 0, and the associated (axial- vector) Nocther current is 

5 H 

5 (Of^(pb) 

1 

= 2 (f^^CTfc - « d'^TTb) dabc (c^c + ilTc)+ h-C 

= dabc{<Tbd^Tic~^bd^a,) . (A2) 

Inserting this into the PCAC relation, 

(OlJrK)EEzp'^/, , (A3) 

where fa is the decay constant corresponding to the field tTq, and shifting the scalar fields by their 
vacuum expectation values, (Ja ^ <^a + ^a, one obtains 

fa = daab &b , (A4) 

where one sums over the index b but not over a. 

In the case that ao, CTg ^0, one obtains for the pion and kaon decay constants 
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/tt = /l — dlla O'a 
Ik ^ fi — diia CTq 



2 

2 _ 1 



CT8 



(A5a) 
(A5b) 



In the case that erg = 0, this shiipHfics to /^r = = -^Z 2 / 3 ctq . 
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